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We study the problem of self-assembly of nanoparticles (NPs) into finite mesoscopic structures 
with a programmed local morphology and complex overall shape. Our proposal is to use NPs 
directionally-functionalized with DNA (dfNP), as building blocks. The combination of directionality 
and selectivity of interactions allows one to avoid any unwanted metastable configurations that 
would commonly lead to slow self-assembly kinetics even in much simpler systems. With numerical 
simulations, we show that a variety of target mesoscopic objects, can be designed and self-assembled 
in near perfect yield. They include cubes, pyramids, boxes and even an Empire State Building model. 
We quantitatively describe the observed kinetics with a simple model, and put forward a compact 
set of design principles for dfNP-based mesostructures. 



PACS numbers: 81.16.Dn,87.14.gk,81.07.-b 

Impressive progress has been made recently in using 
the molecular recognition properties of DNA to direct the 
self-assembly of nanoparticles [TJ [2] into a variety of peri- 
odic nanostructures (3H5] . Many experimental [SHIT] and 
theoretical [T2r[To] studies have aided our understanding 
of these systems and helped in formulating design rules. 
However, the NP-DNA superlattices that were in the fo- 
cus of most of the research to date, represent only a small 
subset of the plausible morphologies for this class of sys- 
tems. One may expect that DNA hybridization would 
enable one to design arbitrary nanoparticle clusters with 
a prescribed overall architecture. In this work, we ex- 
plore a methodology for "programming" such NP-DNA 
mesoscopic structures, and test it with computer simula- 
tions. We demonstrate that the slow kinetics and errors 
typically associated with unwanted metastable states can 
be completely avoided for nanoparticles which are direc- 
tionally functionalized with DNA. 

Conceptually, we are looking for the solution of the 
inverse problem in self-assembly: To identify the set of 
constituent particles that will form the desired structure 
in a robust, error-free manner. In a finite set of colloids 
isotropically functionalized with DNA, one can design the 
pairwise DNA-mediated interactions so that a particular 
target cluster will be its ground state [T7HU]- However, 
it has been shown both theoretically and experimentally 
that the space of metastable configurations in such a sys- 
tem grows rapidly with the cluster size, making the over- 
all kinetics very slow. Furthermore, in order to program a 
structure containing n particles one would typically need 
to use n distinct DNA sequences to encode each indi- 
vidual particle differently, and isolate a complete set of 
the particles prior to the self-assembly. In the theoretical 
paper that precedes this work [2D], we have argued that 
these problems are avoided by the use of directionally- 
functionalized nanoparticles (dfNPs) (see Fig. [IJ). It has 
been argued that the target structure is kinetically pre- 
ferred if the partially built mesostructure is sufficiently 
rigid to avoid undesired bonds. We have also estimated 




FIG. 1. (color online). Nanoparticles isotropically functional- 
ized with complementary single-stranded DNA self-assemble 
into periodic nanostructures like a body centered cubic crys- 
tal (unit cell drawn in black). When the nanoparticles are 
directionally functionalized by attaching DNA strands at spe- 
cific locations, the particles can be designed to self-assemble 
into finite-size mesoscopic architectures like the Empire State 
Building. 

the upper bound of the number of DNA sequences needed 
to program a particular structure. In this work, we use 
computer simulations to move from that abstract the- 
oretical foundation towards more specific and practical 
realizations. 

The dfNPs are more complex than regular DNA- 
functionalized particles due to the fact that single- 
stranded DNA (ssDNA) with given sequences are at- 
tached at well-defined positions rather than randomly. 
Recent experiments indicate that achieving a greater con- 
trol over DNA functionalization and eventually building 
arbitrary dfNPs is within reach. For instance, Suzuki 
et al. 21J have made particles with two distinct ssDNA 



attached at a well-defined separation. Kim et al. [22] 
used electrostatic repulsion to form dfNPs with up to 
6 ssDNA, and colloidal particles with controlled DNA 
surface patches have also been synthesized [53]. Theo- 
retically, we proposed several pathways to self-assemble 
dfNPs jTTJ [20] , inspired by advances in DNA nanotech- 
nology 24-26 . These building blocks, in addition to 
the mutual recognition of the neighboring particles, pro- 
vide control over valence (number of bonds per par- 
ticle) and, to a large extent, bond directionality. In 
other words, dfNP combine the advantages of two ma- 
jor classes of building blocks proposed for programmable 
self-assembly: high selectivity of binding typical for 
isotropic DNA-functionalized particles and the direction- 
ality of interactions similar to that of patchy colloids [27] . 

One can also make a parallel between our approach 
and the booming field of DNA nanotechnology. Simi- 
lar to DNA origami [2S] and most recent DNA Lego [2B] 
approaches, our objective is to create a mesoscopic 3D 
object with predetermined architecture, and similar to 
DNA tiling and algorithmic self-assembly [2U [23 |2T)] . 
we rely on selective directional interactions to program 
the structure. However, in our case nanoparticles of 
an arbitrary type are naturally incorporated in the de- 
sign, which greatly expands the range of potential appli- 
cations. These NP-DNA mesostructures require a new 
set of design rules (discussed below) , and they provide a 
generic solution to the inverse self-assembly problem. 

Our approach will be to engineer a number of in- 
creasingly complex mesostructures, and to use Brown- 
ian dynamics simulations to verify that they are success- 
fully self-assembled from randomly dispersed dfNPs. In 
these simulations, the particle volume fraction is <f> = 
(4/3)ira 3 N/L 3 ~ 10 -3 , where a is the particle radius, N 
in the number of particles and L is the side of the cubic 
simulation box. The equilibrium bond length between 
particles is taken as I = a/2. All times reported here are 
normalized by l/k a C, where k a is the association rate 
constant for nanoparticles functionalized with a single 
DNA strand each and C = N/L 3 . The details of our 
simulation methodology are found in the Supplementary 



Material (SM) 



We begin by considering the simple case of forming a 
finite cluster where the particles are arranged on the ver- 
tices of a cube. This implies eight particles per structure 
and 12 bonds since the particles will be designed such 
that the bonds form along the edges of the cube. The 
most complicated and robust design would be to encode 
each bond with a unique pair of complementary DNA 
strands leading to 12 complementary pairs and 8 parti- 
cle types (Design 1 in Fig. 2(a)). We can also use the 
underlying symmetry of the target structure to reduce 
the number of DNA sequences and particle types. This 
process resulted in 5 different designs with gradually re- 
ducing complexity as shown in Fig. 2(a). 

Simulations have been conducted to study the self- 
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FIG. 2. (color online), (a) Five designs with the NPs arranged 
on the vertices of a cube. For each design the unique comple- 
mentary ssDNA pairs or bonds are drawn with a unique color 
and particles of the same color are identical. Design 1: 12 
complementary pairs (8 particle types), Design 2: 4 comple- 
mentary pairs (4 particle types), Design 3: 3 complementary 
pairs (2 particle types), Design 4: 1 complementary pair (2 
particle types), Design 5: 1 self-complementary pair (1 parti- 
cle type), (b) Final configuration of the system with Design 
3. This simulation was conducted at <f> = 10~ 2 for illustrative 
purposes, (c) Percent yield of cubes versus time for the five 
designs with cj> — 10~ 3 . 



assembly of each of the five cube designs. At initial- 
ization, the particles are assigned random positions and 
orientations. The percent yield of cubes is plotted versus 
time in Fig. 2 for all of the designs. Designs 1-3 form 
in near perfect yield on the time scale of the simulation. 
While Design 3 has only 3 unique complementary pairs it 
leads to nearly the same yield and kinetics as the design 
with 12 pairs (Video 1 of SM 30J). Design 4 is produced 
in a yield of less than 40%, and design 5 leads to the 
formation of tetrahedra and 6-particle clusters in high 
yield with very few cubes. From the experimental point 
of view, Design 3 is optimal because it requires the fewest 
number of unique DNA sequences to produce the target 
cluster quickly and in near perfect yield. We conclude 
that symmetry can be used to reduce the complexity of 
the system, subject to two limitations: (i) neighbors have 
to be of different types and (ii) all bonds for a given 
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FIG. 3. (color online), (a) Pyramids can be formed in a generation-by-generation manner by using a central seed particle. 
While the pyramid is composed of 60 bonds and n = 30 particles, our encoding scheme needs only 15 unique bonds and 9 
particle types. The key to the generation-by-generation design principle can be understood by considering a green particle. 
When an unbound green particle encounters either an unbound red or yellow particle it can form at most one bond which will 
break at some later time. However, when Generation 1 is in place and a green particle bonds to either a red or yellow particle 
it will very quickly form two additional bonds keeping it in place, (b) Final configuration of the pyramid system with excess 
particles and (j> = 0.0025. (c) Percent yield of pyramids versus time with and without excess particles and cj> = 0.0025. Fits to 
Eq. 2 are shown as a dashed black line for the case of extra particles and in orange for the case without. 



particle should be encoded by different sequences. This 
simple rule holds well beyond the simplest cubic clusters, 
and remains valid for most mesostructures studied below. 
In real experiments, a rough estimate of the characteris- 
tic time scale is l/k a C ~ 10 s, which indicates that the 
cubes form in approximately 10 minutes. 

To fabricate structures more complicated than cubes 
it becomes important to employ a strategy that controls 
both the process of nucleation and growth. Ideally, one 
wants to operate in the regime where the number of nu- 
cleation sites or dynamically-growing clusters matches 
the total number of possible clusters, with all other par- 
ticles remaining unbound until joining a cluster. Such 
conditions can very nearly be realized by forming the 
mesostructure in a generation-by-generation manner as 
shown for a pyramid in Fig. 3(a). The idea is to use a 
seed or initiator particle to template the assembly. Par- 
ticles in the first generation form very strong bonds to 
the seed particle. This can be accomplished experimen- 
tally by using DNA sequences with a high CG content. 
Particles in higher generations then add easily because 
they can quickly form multiple bonds with the particles 
that are already in place. This approach suppresses un- 
wanted nucleation which may interfere with a growing 
cluster, leading to mutant structures and slow kinetics. 
The success of the generation-by-generation approach is 



demonstrated in Fig. 3(b) where pyramids are shown to 
self-assemble in perfect yield (see Video 2 of SM 30J). 

The yield of pyramids as a function of time is shown 
in Fig. 3(c) for two cases: when the system composi- 
tion correspond to exactly 4 pyramids, and when extra 
particles (not seed) have been added. The kinetics are 
significantly faster when extra particles are present but 
both cases lead to perfect yield. When the clusters are 
in an intermediate state of growth, a simple description 
of the time evolution of the concentration Cj of particles 
of type i that have yet to join the clusters is 



(1) 



In Eq. [T] c cx t ra corresponds to the extra particles, and 
k a ~ k a is the corresponding association rate constant. 
In the first case, c ox t r a = 0, the concentration of vacancies 
decays slowly as a power law, C{ ~ 1/fca t, while adding 
extra particles leads to sharper exponential relaxation, 
Cj ~ exp(— fc*Cext ra t). The observed yield corresponds to 
the fraction of fully formed clusters, which can be esti- 
mated as 
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FIG. 4. (color online), (a) An orthorhombic box structure self-assembles in a generation-by-generation manner when templated 
by a linear seed. The 114 bonds and n = 63 particles of the structure are encoded using only 11 unique bonds and 13 particle 
types. Note that the length of the structure can be increased with a minimal increase in the number of particles types, (b) The 
configurations at the beginning and end of the low temperature stage are shown for the box structure. Note that four linear 
seeds are visible in the starting configuration, (c) Percent yield versus time for the box structure during the low temperature 
phase, (d) An Empire State Building structure composed of 210 particles (51 types) and 476 bonds (96 unique bonds). 



where is number of particles of type i in a cluster and 
n is the finite cluster size. The data in Fig. 3(c) show 
good agreement with this simple model, indicating that 
the self-assembly is limited primarily by diffusion. 

While symmetry is helpful in reducing the design com- 
plexity, there are many finite-size objects with a repeat- 
ing motif which lack symmetry in a strict mathematical 
sense. For instance, in the 7-layer box structure shown 
in Fig. 4(a), the translational symmetry is broken due 
to its finite size. Does this imply that each of the layers 
has to be built out of distinct set of particles? In the 
proposed design, we avoided this by using a linear seed 
made of 7 distinct particles, that effectively encodes the 
size and ID coordinate along the cluster. The rest of 
the box (56 particles) is built around the seed with only 
6 distinct particle types. This idea is potentially very 
powerful for scalability of our methodology. To illustrate 
this, we have combined the above strategies and design 
rules into constructing an Empire State Building model 
made of 210 dfNP (51 types) as shown in Fig. 4(d) (see 
Video 3 of SM 30 ) . While it is unlikely that such a struc- 
ture will be realized experimentally in the near future, it 



serves as a remarkable demonstration of the potential of 
NP-DNA self-assembly with directional bonding. This 
ability to produce mesostructures with complicated pro- 
grammed architecture, in a parallel and nearly error-free 
manner, can eventually lead to a great variety of practi- 
cally important applications in the fields of metamateri- 
als, nanoplasmonics 31 a and nanomedicinc. 
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